arXiv:1508.06377v2 [cs.SY] 21 Mar 2016 


Performance Analysis and Coherent Guaranteed Cost Control for 
Uncertain Quantum Systems Using Small Gain and Popov Methods 

Chengdi Xiang, Ian R. Petersen and Daoyi Dong 


Abstract —This paper extends applications of the quantum 
small gain and Popov methods from existing results on robust 
stability to performance analysis results for a class of uncertain 
quantum systems. This class of systems involves a nominal lin¬ 
ear quantum system and is subject to quadratic perturbations in 
the system Hamiltonian. Based on these two methods, coherent 
guaranteed cost controllers are designed for a given quantum 
system to achieve improved control performance. An illustrative 
example also shows that the quantum Popov approach can 
obtain less conservative results than the quantum small gain 
approach for the same uncertain quantum system. 

I. INTRODUCTION 

Due to recent advances in quantum and nano-technology, 
there has been a considerable attention focusing on research 
in the area of quantum feedback control systems; e.g., [1]- 
[24], In particular, robust control has been recognized as 
a critical issue in quantum control systems, since many 
quantum systems are unavoidably subject to disturbances and 
uncertainties in practical applications; e.g., [1], [2] and [3]. 
A majority of papers in the area of quantum feedback control 
only consider the case in which the controller is a classical 
system. In this case, analog and digital electronic devices 
may be involved and quantum measurements are required 
in the feedback loop. Due to the limitations imposed by 
quantum mechanics on the use of quantum measurement, re¬ 
cent research has considered the design of coherent quantum 
controllers to achieve improved performance. In this case, the 
controller itself is a quantum system; e.g., [1], [4] and [5]. 
In the linear case, the quantum system is often described 
by linear quantum stochastic differential equations (QSDEs) 
that require physical realizability conditions in terms of 
constraints on the system matrices to represent physically 
meaningful systems. 

As opposed to using QSDEs, the papers [6], [7] have 
introduced a set of parameterizations ( S , L, //) to represent 
a class of open quantum systems, where S' is a scattering 
matrix, L is a vector of coupling operators and H is a 
Hamiltonian operator. The matrix S, together with the vector 
L , describes the interface between the system and the field, 
and the operator H represents the energy of the system. The 
advantage of using a triple (S, L, H) is that this framework 
automatically represents a physically realizable quantum 
system. Therefore, in this paper, a coherent guaranteed cost 
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controller is designed based on ( S,L,H ) and the physical 
realizability condition does not need to be considered. 

The small gain theorem and the Popov approach are two 
of the most important methods for the analysis of robust 
stability in classical control. The paper [8] has applied the 
small gain method to obtain a robust stability result for 
uncertain quantum systems. This result gives a sufficient 
condition for robust stability in terms of a strict bounded real 
condition. The small gain method has also been extended 
to the robust stability analysis of quantum systems with 
different perturbations and applications (e.g., see [9], [10], 
[11] and [12]). The paper [13] has introduced a quantum 
version of the Popov stability criterion in terms of a fre¬ 
quency domain condition which is of the same form as the 
classical Popov stability condition (e.g., see [25]). The Popov 
approach has also been used to analyze the robust stability 
of an optical cavity containing a saturated Kerr medium 
[14], Also, the paper [13] has shown that the frequency 
domain condition obtained in [13] is less conservative than 
the stability result using the small gain theorem [8]. 

In this paper, we extend the quantum small gain method 
in [8] and the Popov type approach in [13] from robust 
stability analysis to robust performance analysis for uncertain 
quantum systems. We assume that the system Hamiltonian 
H can be decomposed in terms of a known nominal Hamil¬ 
tonian Hi and a perturbation Hamiltonian H 2 , i.e., // = 
H 1 + H 2 - The perturbation Hamiltonian // 2 is contained in 
a set of Hamiltonians W. We consider uncertain quantum 
systems where the nominal system is a linear system and 
the perturbation Hamiltonian is quadratic. Moreover, a co¬ 
herent controller is designed using the small gain approach 
and the Popov approach for the uncertain quantum system, 
where a guaranteed bound on a cost function is derived in 
terms of linear matrix inequality (LMI) conditions. Although 
preliminary versions of the results in this paper have been 
presented in the conference papers [22] and [23], this paper 
presents complete proofs of the main results and modifies the 
example in [22] for a consistent performance comparison of 
the proposed two methods. 

The remainder of the paper proceeds as follows. In Sec¬ 
tion [II] we define the general class of quantum systems 
under consideration and specify the nominal system as a 
linear quantum system. We then present a set of quadratic 
Hamiltonian perturbations in Section III In Section |IV[ a 
performance cost function for the given system is defined. 
Moreover, a small gain approach and a Popov type method 
are used to analyze the performance of the given system. 
In Section [V] a quantum controller is added to the original 



system to stabilize the system and also to achieve improved 
performance. Also, corresponding approaches are used to 
construct a coherent guaranteed cost controller in terms of 
LMI conditions. In Section VI an illustrative example is 
provided to demonstrate the method that is developed in 
this paper. A performance comparison between these two 
methods is also shown in the illustrative example. We present 
some conclusions in Section IVII1 


II. QUANTUM SYSTEMS 


In this section, we describe the general class of quantum 
systems under consideration, which is defined by parameters 
( S,L,H ). Here H = Hi + H 2 , f?i is a known self-adjoint 
operator on the underlying Hilbert space referred to as the 
nominal Hamiltonian and H 2 is a self-adjoint operator on 
the underlying Hilbert space referred to as the perturbation 
Hamiltonian contained in a specified set of Hamiltonians 
W; e.g., [6], [7], The set W can correspond to a set of 
exosystems (see, [6]). The corresponding generator for this 
class of quantum systems is given by 


g(x) = -i[x,H\ + c(x) (i) 

where C(X) = \tf[X,L\ + \[L\X]L. Here, [X,H] = 
XH — HX describes the commutator between two operators 
and the notation t refers to the adjoint transpose of a vector 
of operators. Based on a quantum stochastic differential 
equation, the triple ( S,L,H ), together with the correspond¬ 
ing generators, defines the Heisenberg evolution X(t) of an 
operator X [6]. The results presented in this paper will build 
on the following results from [13]. 

Lemma 1: [13] Consider an open quantum system de¬ 

fined by (S, L. H) and suppose there exist non-negative self- 
adjoint operators V and W on the underlying Hilbert space 
such that 

G(V) + W < A (2) 


where A is a real number. Then for any plant state, we have 
lim sup ^ f {W(t))dt< X. (3) 

T-t 00 T J 0 

Here W(t ) denotes the Heisenberg evolution of the operator 
W and (•) denotes quantum expectation; e.g., see [13] and 
[ 6 ]. 

In this paper, the nominal system is considered to be 
a linear quantum system. We assume that Hi is in the 
following form 


Hi 


\[ «' « T ]M 



(4) 


where M £ C 2nx2n is a Hermitian matrix and has the 
following form with Mi = M\ and M 2 = M£ 


M = 


Mi M2 

M* M* 


(5) 


Here a is a vector of annihilation operators on the underlying 
Hilbert space and a # is the corresponding vector of creation 
operators. In the case of matrices, the notation '' refers to 


the complex conjugate transpose of a matrix. In the case 
of vectors of operators, the notation # refers to the vector 
of adjoint operators and in the case of complex matrices, 
this notation refers to the complex conjugate matrix. The 
canonical commutation relations between annihilation and 
creation operators are described in the following way 


' 

a 


a 

t ■ 


a 
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. a * . 


. a * . 
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—1 

<3 
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<3 
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( 6 ) 


where J = 


I 0 
0 -I 


[ 2 ]. 


The coupling vector L is assumed to be of the form 


L = [ Ni N 2 


a 

a# 


= N 


7 # 


(7) 


where N\ £ C mxrl and N 2 £ C mxn . We also write 


L 

= N 

a 


' Ni 

n 2 

a 

L* 

. a * . 


. N* 

N f _ 

. a *. 


When the nominal Hamiltonian H is a quadratic function 
of the creation and annihilation operators as shown in (|4]) 
and the coupling operator vector is a linear function of 
the creation and annihilation operators, the nominal system 
corresponds to a linear quantum system (see, [2]). 

We consider self-adjoint “Lyapunov” operators V of the 
form 


V = [ at a T ] P 



(9) 


where P £ C 2nx2n is a positive definite Hermitian matrix 
of the form 


P = 


Pi P 2 

p# p# 

i 2 ‘\ 


( 10 ) 


We then consider a set of non-negative self-adjoint oper¬ 
ators V defined as 


V = 


V of the form such that P > 0 is a 


Hermitian matrix of the form (101 


(ID 


III. PERTURBATIONS OF THE HAMILTONIAN 


In Section [Tl| we introduced the nominal linear quan¬ 
tum system. This section defines the quadratic Hamiltonian 
perturbations (e.g., see [8], [18]) for the quantum system 
under consideration. We first define two general sets of 
Hamiltonians in terms of a commutator decomposition, and 
then present two specific sets of quadratic Hamiltonian 
perturbations. 








































A. Commutator Decomposition 

For the set of non-negative self-adjoint operators V and 
given real parameters 7 > 0 ,5 > 0 , a particular set of 
perturbation Hamiltonians Wi is defined in terms of the 
commutator decomposition 

[V,H 2 ] = \V,zi]w-wi[z,V\ (12) 


for V £ V, where w and 2 are given vectors of operators. 
W 1 is then defined in terms of sector bound condition: 


v^w < —77 z 4- 5. 

T 


We define 

Wi = 


H 2 : 3 w, 2 such that (13) is satisfied 
and ( 12 ) is satisfiedv V £ V 


(13) 


(14) 


B. Alternative Commutator Decomposition 


Given a set of non-negative operators V, a self-adjoint 
operator Hi, a coupling operator L, real parameters /3 > 0 
7 > 0, and a set of Popov scaling parameters 0 C [0, 00 ), 
we define a set of perturbation Hamiltonians VV ’2 in terms of 
the commutator decompositions [13] 


[V-9H 1 ,H 2 ] = [V - 0H u z^]w - w^[z,V - eH x ], 
£(H 2 ) < C(z^)w + w^C(z) + j3[z, L]1 [z, L) 


(15) 


for V £ V and 9 £ 0, where w and 3 are given vectors of 
operators. Note that © and ( fl5j ) correspond to a general 
quadratic perturbation of the Hamiltonian. This set W 2 is 
then defined in terms of the sector bound condition 

1 


(w - z)l('U) - z) 

7 7 


< 


-z'z. 


7 


We define 


W 2 = 


H 2 > 0 : 3 w,z such that (15) and (16) 


are satisfied \/ V £ V~U £ 0 


(16) 


(17) 


C. Quadratic Hamiltonian Perturbation 

We consider a set of quadratic perturbation Hamiltonians 
that is in the form 


Hi = \[ c e ]a 


C 

c # 


(18) 


where £ = E\a + E 2 a# and A £ C 2mx2m is a Hermitian 
matrix of the form 


A = 


At A 2 
A* Af 


(19) 


with Ai = A{ and A 2 = A 2 . 

Since the nominal system is linear, we use the relationship: 


2 = 


Then 


c 

El 

e 2 

a 

C # J “ [ 

Ef 

E* 

. a * 


a f 

a T ] 

E^AE 


= E 



a* ■ 
( 20 ) 


( 21 ) 


When the matrix A is subject to the norm bound 

||A|| < ( 22 ) 

7 

where II. 11 refers to the matrix induced norm, we define 


W 3 = 


H 2 of the form ( p~8] > and ( | 1 9[ > such that 


condition 


is satisfied 


(23) 

In [ 8 ], it has been proven that for any set of self-adjoint 
operators V, 

W 3 CW 1 . (24) 


When the matrix A is subject to the bounds 


0 < A < -I, 
7 


(25) 


we define 


W 4 = 


J II 2 of the form ( [T 8 | and ( fl9| such that 


condition 


is satisfied 


(26) 

In [13], it has been proven that if [ 2 , L] is a constant vector, 
then for any set of self-adjoint operators V, 

W 4 C W 2 . (27) 

IV. PERFORMANCE ANALYSIS 

In this section, we provide several results on performance 
analysis for quantum systems subject to a quadratic per¬ 
turbation Hamiltonian. Also, the associated cost function is 
defined in the following way: 

rT 


J = lim sup 

T—too J- 


([ a 1 a T ] R 


7# 


where R > 0. We denote 


W = [ 


it 


R 


1 # 


)dt (28) 


(29) 


In order to prove the following theorems on performance 
analysis, we require some algebraic identities. 

Lemma 2: (See Lemma 4 of [ 8 ]) Suppose V £V,H\ is 
of the form Q and L is of the form ([7]). Then 

it 


\V,Hi] = 


7 # 


(PJM - MJP) 


7# 


(30) 


£ ^ = -2 


it 


7# 


Tr(P JTVl 


(N^JNJP + PJN^JN) 
NJ ), 


I 0 
0 0 


a 


Lemma 3: For V £ V and 2 defined in (20]), 


= 2 JP 


7# 


[ 2 , V] = 2 EJP 


7# 


7# 


(31) 

(32) 

(33) 


[V,z*][z,V]= 4 


it 


7# 


PJE^EJP 


7# 


(34) 

















































z* z = 


it 


T# 


E ] E 


Proof: The proof follows from Lemma [2] 


(35) 

□ 


Lemma 4: (See Lemma 5 of [13]) For 2 defined in ( |20| ) 
and L being of the form ((TJ), 


[z,L\ = [E 

is a constant vector, where 

£ = 


a 


a 

. a * . 

,N 

. a * . 


] = EJT,N t 


(36) 


0 / 
/ 0 

Lemma 5: (See Lemma 6 of [13 


(37) 


) For 2 defined in 


Hi defined in (|4ji and L being of the form (|7j, we have 


- i[z, Pi] + £{z) = E(—iJM - -JN*JN) 


and 


i[z,V} = 2iEJP 


t# 


(38) 


(39) 


Now we present two theorems (Theorem 1 and Theorem 2) 
which can be used to analyze the performance of the given 
quantum systems using a quantum small gain method and a 
Popov type approach, respectively. 

A. Performance Analysis Using the Small Gain Approach 

Theorem 1: Consider an uncertain quantum system 
(S, L, H ), where H = Hi + H 2 , Hi is in the form of (pj, L 
is of the form ([7J) and P 2 £ W 3 . If F = —iJM— 3 JJAPJN 
is Hurwitz, 

E ] E 


Ftp + PF + + R 2PJFt 

2 EJP -I/t 2 


< 0 (40) 

has a solution P > 0 in the form of ( [T()| and t > 0 , then 

J = lim sup f ( W(t))dt 
T —>00 -I Jo 


= lim sup — ([ at a T ] R 


t—> c 




)dt < A H—— 


where 


A = Tr(PJiVt 


I 0 
0 0 


NJ). 


T 

(41) 


(42) 


In order to prove this theorem, we need the following lemma. 


Lemma 6 : Consider an open quantum system (.S', L. H) 
where H = Hi + H 2 and H 2 £ Wi, and the set of non¬ 
negative self-adjoint operators V. If there exists a V £ V 
and real constants A > 0, r > 0 such that 

- i[V, Pi] + C(V) + r 2 [F, 2 t][ 2 , V] + + W < A, 

(43) 

then 

1 r'l' s 

limsup— / (W(t))dt < A 4 —> 0. (44) 

T-voo 7 J Q t- 


Proof: Since V £ V and Il 2 £ VV’i, 

g(V) = -i[V,H 1 ]+/:{V)-i[V,z^w + iw^[z,V]. (45) 


Also, 

0 < (r[V, 2 t]- u; t)(r[V, 2 ;t] -tt)t)t 

T ^ 


l 

T 


= t 2 [V, z^] [z , V] + i\V, z^]w — iw^[z, V] + 


w'w 


(46) 


Substituting © into © and using the sector bound 
condition ( fl3] l, the following inequality is obtained: 

1 


G(V) < -z{V,Hi}+C(V)+t 2 [V,z^}[z : V} + 
Hence, 


„ n Z Z~\ . 

7 Z T Z T z 


g(v) + w<\ + 


-2 ‘ 


S_ 

(47) 

(48) 


Consequently, the conclusion in the lemma follows from 
Lemma [l] □ 

Proof of Theorem [7] Using the Schur complement [26], 
the inequality (|40|) is equivalent to 


F f P + PF + 4t 2 PJE^ EJP + 


FtF 

7 2 t 2 


+ R< 0. (49) 


If the Riccati inequality ( [49] ) has a solution P > 0 of the 
form (TO] and r > 0, according to Lemma [2] and Lemma [3] 
we have 


- t[v, Pi] + C(V) + t 2 [V, 2t] [2, V] + ^zh + 


7 2 r 2 


W = 


a 

f / F^P + PF + 4t 2 PJE^EJP \ 

a 

. a * . 


. a * . 


Tr(PJ7V t 


/ 0 
0 0 


NJ). 


(50) 

Therefore, it follows from ( |40| that condition ( |43] > is satisfied 
with 

' I 0 
0 0 


A = Tr {PJN^ 


NJ) > 0. 


(51) 


Then, according to the relationship (24 1 and Lemma [ 6 ] 
have 


lim sup — 
T—foo T 


(W(t))dt 


= lim sup — ([ a) a T ] R 

T^>oo J Jo 




)dt < X -\—. 


(52) 

□ 


B. Performance Analysis Using the Popov Approach 

Theorem 2: Consider an uncertain quantum system 
(S, L , P), where H = Hi + P 2 , Hi is in the form of (pji, L 
is of the form (j7j) and il 2 £ Wi. If F = —iJM — fjhPJN 
is Hurwitz, and 


PF + Ftp + P -2fPJFt + pt + 0FtFt 
2iEJP + E + dEF - 7 1 


< 0 
(53) 








































has a solution P > 0 in the form of ( fl0| ) for some 9 > 0, 
then 


J = lim sup f ( W(t))dt 

T—t oo J- JO 

i r T 

= lim sup — I ([ a/ a T ] R 

T —700 T J q 

where 

A = Tr(PJ7V t 


7# 


)dt < A (54) 


I 0 
0 0 


NJ) 


49 


N*YJE^EJYN T . 


7 


Then 


1 ( T 

lim sup— / (W(t))dt < A. (57) 

T—7oo 7 JO 


Here W(i) denotes the Heisenberg evolution of the operator 
W and (•) denotes quantum expectation. 

Proof of Theorem [Aj- Using the Schur complement, (531 is 
equivalent to 

PF + F^P + 1 (-2 iPJE* + F* + GF^EO 

7 (58) 

x (2 iEJP + E + 9EF ) + F < 0. 

According to Lemma [2] and Lemma [5] we have 

- i[V, Hi] + C(V) + ~(i[z, V - 9Hi] + 9C(z) + *)* 

7 

x (i[z, V-9H i] + 9C{z) + z) + — [z, L]t [z, L] + W 

7 

PF + Ftp 

+ i(-2iPJFt + Ft + 6»ptFt) 

x (2 iEJP + E + 9EF) + P 


7# 


t 


# 


+ Tr(P JTVt 


I 0 
0 0 


49 


NJ) + —N*Y 1 JE^EJT,N t . 


(59) 


From this and using the relationship ( |27j ). Lemma |4] and 
Lemma [7] we obtain 


i r T 

limsup— / ( W(t))dt 

T->oo 7 Jq 


= lim sup — 

T —foo 7 1 , 


nt 


7# 


7? 


7# 


)dt 


(60) 


where 

A = Tr(PJAtt 


I 0 
0 0 


49 


NJ) + —7V # SJF t FJ£iV T . 


(61) 

□ 


(55) 

In order to prove this theorem, we need the following lemma. 

Lemma 7: (See Theorem 1 of [13]) Consider a set of non¬ 
negative self-adjoint operators V, an open quantum system 
(S. L. H) and an observable W, where // = H \ + II 2 and 
H 2 € W 2 defined in ( fl7| ). Suppose there exists a V € V and 
real constants 9 £ 0, A > 0 such that 

- i[v, Pi] + ay) + ~{i[z, v - 9 h 1 ] + ec{ z ) + z) t 
7 

x (i[z, V — 9Hi] + 9C{z) + z) + 9/3[z, L]^[z, L\ + W < A. 

(56) 


V. COHERENT GUARANTEED COST 
CONTROLLER DESIGN 

In this section, we design a coherent guaranteed cost 
controller for the uncertain quantum system subject to a 
quadratic perturbation Hamiltonian to make the control sys¬ 
tem not only stable but also to achieve an adequate level of 
performance. The coherent controller is realized by adding a 
controller Hamiltonian H 3 . // . is assumed to be in the form 


H,=\{ a 


K 


7# 


(62) 


where K £ C 


2nx2n 


is a Hermitian matrix of the form 


K = 


Ki I < 2 
K* Kf 


(63) 


and K-\ = I\\, K 2 = K) . Associated with this system is 
the cost function J 


J = limsup — 

T—>00 J- , 


it 


7 # 


(R + pK 2 ) 


7# 


where p £ (0, 00 ) is a weighting factor. We let 
it 


W = 


7# 


(R + pK 2 ) 


7# 


)dt (64) 


(65) 


The following theorems (Theorem 3 and Theorem 4) 
present our main results on coherent guaranteed cost con¬ 
troller design for the given quantum system using a quantum 
small gain method and a Popov type approach, respectively. 

A. Coherent Controller Design Using the Small Gain Ap¬ 
proach 

Theorem 3: Consider an uncertain quantum system 
(S,L,H), where H = H\ + H 2 + // 3 , Hi is in the form 
of ([4]), L is of the form 0. II > £ VV .3 and the controller 
Hamiltonian H 3 is in the form of (62). With Q = P~ l , 
Y = KQ and F = -iJM - \JNUN, if there exist a 
matrix Q = q* I (q is a constant scalar and / is the identity 
matrix), a Hermitian matrix Y and a constant r > 0, such 
that 


<0 ( 66 ) 


where A = qF t + Fq + iYJ — iJY, then the associated cost 
function satisfies the bound 

rT 


A + 4 t 2 JE^EJ 

Y 

gf?3 

qE* 

Y 

- I/P 

0 

0 

qRi 

0 

-I 

0 

qE 

0 

0 

-7 2 t 2 I _ 


limsup - / 
T—700 1 Jo 

= limsup- 
T—f 00 * 

r (5 


(' W(t))dt 


it 


7# 


(R + pK 2 ) 


7# 


)dt m 






































where 


A = Tr {PJN^ 


I 0 
0 0 


NJ). 


( 68 ) 


Proof: Suppose the conditions of the theorem are satisfied. 
Using the Schur complement, ( | 66 | is equivalent to 

A + At 1 JEfEJ + pYY qRi gPt 

qRl -I 0 < 0 . 

qE 0 -7 2 t 2 I 

(69) 


<0 (70) 


Applying the Schur complement again, it follows that (69i 
is equivalent to 

A + 4 r 2 JE^EJ + pYY + q 2 R gPt 

qE -7 2 t 2 I 

and ( |70| > is equivalent to 

qF f + Fq + iYJ - iJY + 4 t 2 JE^EJ 

,E^E 


+ pYY + q 2 (EUA +R ) <0. 


(71) 


'7 2 t 2 


Substituting Y = Kq = qK t into (71 1 , we obtain 

q(F - iJK )t + {F- iJK)q + 4 t 2 JE^EJ 

2 .ptp n 2 „ (72) 

+ 9 (— 7 ,— 7 + P -T pK ) < 0. 

7~' 

Since P = Q _1 , premultiplying and postmultiplying this 
inequality by the matrix P, we have 

(P - zJP) t P + P(F - iJK) + 4 t 2 PJE^EJP 

ptp „ 2 n ( 73 ) 

H—^^ + P + pK < 0. 

It follows straightforwardly from (731 that F — iJK is 
Hurwitz. We also know that 

- i[V, H x + H 3 ] + C(V) + t 2 [V, z*][z, U] + -^zh + W 

(P - UK^P + P(P - iJK) 
+4 t 2 PJE^EJP + E^E/(j 2 t 2 ) 

+R + pK 2 

I 0 


Tr(PJA rt 


0 0 


NJ). 


According to the relationship ([24]) and Lemma [ 6 ] we have 


1 r 

limsup — / ( W(t))dt 

T —>00 T J 0 


= limsup- 

T—f 00 . 

^ \ $ 

< A H —~ 


0 

2 roo 


0 


7 # 


(P + pK 2 ) 


r# 


)dt 


(75) 


where 


We know that P = Q 1 = q 1 I and apply the Schur 
complement to inequality (f77[, so that we have 


■£ + 72 
P5 


P3 

-9 


< 0 


(78) 


I 0 
0 0 


where B = Tr(JTVt 
complement again, it is clear that (|78]> is equivalent to 


NJ). Applying the Schur 


B 2 

Si -t 2 0 

P 5 0 -q J 


< 0 . 


Hence, we minimize £ subject to ( |79[ and 
[3] This is a standard LMI problem. 


(79) 


in Theorem 


B. Coherent Controller Design Using the Popov Approach 

Theorem 4: Consider an uncertain quantum system 
(S', L. II ), where II = U\ + II 2 + Il:>,, Hi is in the form 
of <[4]), L is of the form <[T[), H 2 £ W 4 , the controller 
Hamiltonian /7 :i is in the form of (|62|). With Q = P _1 


Y = KQ and P = —iJM — \ JNUN , if there exist a 
matrix Q = q * I (q is a constant scalar and / is the identity 
matrix), a Hermitian matrix Y and a constant 6 > 0, such 
that 


A 

B t 

Y 

qRi ' 

B 

- 7 / 

0 

0 

Y 

0 

-I/p 

0 

_ qRi 

0 

0 

-I 


< 0 


(80) 


where A = Fq + qF t — iJY + iYJ and B = 2iP J + Eq + 
6 EFq — idEJY, then the associated cost function satisfies 



the bound 

a 

1 

a * _ 

lim sup — 


T— >00 T 

(74) 

= lim sup 

T—¥ OO 


. .„ h 

where 

A = Tr(PJ7V t 


I 0 
0 0 


(.R + pK 2 ) 


40 


7# 


)dt < A 

(81) 


NJ) + —N*YJE^EJYN T . 

7 

(82) 

Proof: Suppose the conditions of the theorem are satisfied. 


Using the Schur complement, (80i is equivalent to 


A = Tr(PJN^ 


'I O' 

NJ). 

(76) 

r A + ^-B^B 

Y 

qR^ ' 

0 0 

7 

Y 

-I/p 

0 



□ 

gP3 

0 

-I 


< 0 . 


(83) 


Remark 1: In order to design a coherent controller which 
minimizes the cost bound © in Theorem 3, we need to 


Tr(PJ7V t 


We then apply the Schur complement to inequality ( [83] ) and 
obtain 


< 0 . 


I 0 " 

6 

A+ IB^B + pYY qR^ 

0 0 

NJ) H—^ < e (77) 

l 

toll— 1 

1 


(84) 
















































Also, (84 1 is equivalent to 


Fq + qF t - iJY + iYJ+ 

-(-2 iJE* + qE f + 9qF ] E^ + iOYJE *) 
7 

x (2 iEJ + Eq + OEFq - iOEJY) 

+ q 2 R + pYY < 0. 


Substituting Y 


Kq = qK t into (851, we obtain 


(F - iJK)q + q(F - * JK )* 

+ -(-2 iJE* + qE t + 9q{F - UK^E*) 

7 

x (2 iEJ + Eq + 9E(F - iJK)q) 

+ q 2 (R + pK 2 ) < 0. 


(85) 


( 86 ) 


Since P = Q 1 , premultiplying and postmultiplying this 
inequality by the matrix P, we have 


P{F - iJK) + (F - iJK^P 

+ -(-2 iPJE^ +E f + 9{F - iJK^E^) (87) 

7 

x (2 iEJP + E + 9E(F - UK)) + R + pK 2 < 0. 

It follows straightforwardly from © that F — iJK is 
Hurwitz. We also know that 


-i[v,fli + fr3] + £(v) 

+ ~(i[z, V - 9{H X + H 3 )} + 9C(z) + *) f 
7 

x (i[z, V - 9{H 3 + H 3 )] + 9C{z) + z) 


4 9 


[z,L]'[z,L] + W 




M 


+ Tr {PJN^ 


a# 

I 0 
0 0 


4 9 


NJ) + —N*YJE^EJYN t 


( 88 ) 


where 

M =P(F - UK) + (F — iJK) f P 

+ 2iPJE t + + 

7 \ ) ( 89 ) 

x (2iAJP + E + 9E(F - i JK )) 

+ i? + pK 2 . 


According to the relationship (27 i and Lemma [7] we have 


i r 

limsup— / (W(t))dt 
T—too 7 Jo 


0 

= lim sup — 

T—>■ oo T , 

< A 


it 


i# 


(i? + pit: 2 ) 


T# 


)df 


(90) 


where 

A = Tr(PJA rt 


/ 0 
0 0 


7 VJ) + —N*YJE^EJYN t . 

7 

( 91 ) 


□ 

Remark 2: For each fixed value of 9, the problem is an 
LMI problem. Then, we can iterate on 9 £ [0, oo) and choose 
the value which minimizes the cost bound ( |82| i in Theorem 

El 


VI. ILLUSTRATIVE EXAMPLE 


In order to illustrate our methods and compare their 
performance, we use the same quantum system considered in 
[23] as an example. The system corresponds to a degenerate 
parametric amplifier and its (S, L. H) description has the 
following form 


H=-i{(aJ) 2 -a 2 ),S = I 1 L = 

We let the perturbation Hamiltonian be 


' na. 


H 2 = ^ [ aJ a 1 


1 0.5z 

—0.57 1 



a 


— 1 

e 

_ 1 


(92) 


(93) 


and the nominal Hamiltonian be 


H t = - [ at a 2 
2 L 


-1 0.5i 

-0.5 i -1 



a 


. a * . 


(94) 


so that iii + ii 2 = //. The corresponding parameters con¬ 
sidered in Theorem 1, Theorem 2, Theorem [3] and Theorem 
[4] are as follows: 


M = 
F = 


-1 

0.5*' 

> 7 V = 

\/« 

0 

-0.5* 

-1 

0 

s/n 


-5 +» 
6.5 


0.5 


.E = I 


and 


A = 


1 0.5* 

-0.5* 1 


(95) 


(96) 


To illustrate Theorem 1 and Theorem 3, we consider 
H 2 € W 3 - Hence, 7 = 1 is chosen to satisfy ( |22] ). The 
performance using the small gain approach for the uncertain 
quantum system is shown in Figure 1. In Figure 1, the dashed 
line represents the cost bound for the linear quantum system 
considered in Theorem 1 as a function of the parameter 
n. The solid line shows the system performance with the 
coherent controller designed in Theorem 3. Compared to 
the performance without a controller, the coherent controller 
can guarantee that the system is stable for a larger range 
of the damping parameter k and gives the system improved 
performance. 

Now we illustrate one approach to realizing the desired 
controller. For instance, when k = 4.5, by using the con¬ 
troller design method in Theorem 3, we have the desired 
controller Hamiltonian as 


*. = )[ 


0 

0.5* 



a 


. a * . 


(97) 


This controller Hamiltonian can be realized by connecting 
the degenerate parametric amplifier with a static squeezer as 
shown in Figure 2. This static squeezer is a static Bogoliubov 
component which corresponds to the Bogoliubov transforma¬ 
tion [27], [28]. Also, we have the following definition: 




































Definition 1: (see [27], [28]) 
ponent is a component that ii 

transformation: 


dy(t ) 

— T2 

du(t) 

. dy*{t) 

— jD 

du*{t) 


B\ f?2 


A static Bogoliubov com- 


where B = 


l b r b* j 


,BtJB = J. 


To realize H : >, in (97!, we let the matrix B = 


which satisfies the Bogoliubov condition B^JB = J, and 
k = 77. Therefore, the overall Hamiltonian of the closed 
loop system is H = II\ + // 2 + //.s, which achieves the 
controller design goal. Detailed procedure regarding how to 
get the matrix B can be found in appendix. 



Fig. 1. Guaranteed cost bounds for an uncertain quantum system with a 
controller (solid line) and without a controller (dashed line) using the small 
gain approach 
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Fig. 3. Performance versus 6 


In Figure 4, the dashed line shows the performance for 
the given system considered in Theorem 2 and the solid line 
describes the cost bound for the linear quantum system with 
the coherent controller considered in Theorem 4. As can be 
seen in Figure 4, the system with a controller has better 
performance than the case without a controller. 



Fig. 2. Degenerate parametric amplifier coupled to a static squeezer. 

To illustrate Theorem 2 and Theorem 4, we consider H 2 € 
W 4 . Hence, 7 = 2 is chosen to satisfy ( |25| >. The results 
using the Popov approach are shown in Figure 3 and Figure 
4. Figure 3 demonstrates how to choose the value of 6 . We 
consider the same example as above with k = 3.8 and iterate 
on 9 £ [0,1]. Figure 3 shows the cost bound for this quantum 
system obtained in Theorem 4 as a function of the parameter 
0. It is clear that the minimal cost bound is achieved when 
6 = 0.1. Therefore, we choose 6 = 0.1 for k = 3.8 and use 
a similar method to choose 6 for other values of k. 


Fig. 4. Guaranteed cost bounds for the uncertain quantum system with a 
controller (solid line) and without a controller (dashed line) using the Popov 
approach 

Also, we can observe that the method in Theorem 3 can 
only make the quantum system stable for k > 4 in the 
example. Therefore, compared with the results in Figure 1, 
the Popov method obtains a lower cost bound and a larger 
range of robust stability as shown in Figure 4. This is as 
expected, since the Popov approach allows for a more general 
class of Lyapunov functions than the small gain approach. 

































VII. CONCLUSION 

In this paper, the small gain method and the Popov 
approach, respectively, are used to analyze the performance 
of an uncertain linear quantum system subject to a quadratic 
perturbation in the system Hamiltonian. Then, we add a 
coherent controller to make the given system not only stable 
but also to achieve improved performance. By an illustrative 
example, we have also shown that the Popov method adds 
to a considerable improvement over the small gain method 
in terms of system performance. Future work will include 
the extension of these approaches to nonlinear uncertain 
quantum systems [24]. 


Substituting ([99]) into (lOOi, we obtain that 


B~ 


dA 2 (t) 


VK 

0 

a{t) 

_ dA*(t) _ 


0 

y/k 

a*(t) _ 


dt 


+ 


dA 2 (t ) 
dA*(t) 


Hence, 

(LT 1 - I) 


dA 2 (t) 


\fk 

0 

a(f) 

. dA*(t) _ 


0 


. « # w. 


dt. 

( 101 ) 

We now assume inverse of B~ x — I exists. It follows ( 101 1 
that we can write 


Appendix 

The detailed procedure regarding how to realize the de¬ 
sired controller is shown below. We consider a degenerate 
parametric amplifier (DPA) as an example. Based on the 


(S', L. H) description in (92 1 , we can calculate the following 
quantum stochastic differential equations [1] describing the 
DPA: 


dt 


da(t) 


K 

2 

1 

a(f) 

da*(t ) 


1 

K 

2 . 



v/k 0 

dAi{t) 

0 y/R, 

dA*(t ) 


da(t ) 


K-\-k 

2 

l 

o,(t) 

da*(t) 


l 

K + K, 

2 



0 


dA\{t) 


dA 2 {t) 
dA*(t) _ 


= ( B - 1 - I) 


\-l 


y/k 

0 


0 

y/k 



a{t) 
0 1 ) 


T# 


dt.. 

( 102 ) 


Substituting ( | 102 i into first equation in 

da(t) 
da*{t) 


dA/\t) =y/ka(t)dt + dA\{t). 

We have known that when k = 4.5 and using the controller 
design method in Theorem 3, we have the desired controller 
Hamiltonian as in m 

Next, we show how to realize this controller Hamiltonian 
by connecting this degenerate parametric amplifier with 
a static squeezer as shown in Fig. 2. The corresponding 
quantum stochastic differential equations for this DPA is as 
follows: 


dt 


K-\-k 

2 

1 


a{t) 

l 

K-\-K 

2 


. a* {t) 

y/R 

0 


dAi(t) 

0 

y/k 

_ dAf(t) _ 

y/k 

0 

0 

y/k 

( B 

- 1 - I)~ 

y/k 

0 

‘ 

a(f) 

0 

y/k 

L J 

K-\-k 

2 

1 

' 

a(f) 

l 

K-\-K 

2 


. a# it) 

y/R 

0 


dAi(t) 

0 

y/k 

_ dA*{t) _ 


, we get 
dt 


dt 


-(IT 1 -/)- 1 * J ( > dt. 

0 k J [ a r {t) 

(103) 

According to Definition 1, the B matrix satisfies the 
following relation: 


jbUb = 


I 0 
0 -I 



■ B\ Bl ‘ 


O 


B l B 2 


.4 B ? . 


i 

o 


B* B* _ 


y/k 0 

dA 2 {t) 

(98) 

B\ B/ 

Bi B 2 

0 y/k 

dA*(t ) . 

5 

. ~B\ -Bl _ 

-B* -B* _ 


dA° ut (t) = y/ka(t)dt + dAi(t); 
dA°/ l [t) = y/ka(t)dt + dA 2 [t). 

We have known that this static squeezer is a static Bogoli- 
ubov component which satisfies Definition 1. According to 
the Definition 1, we have the following relation: 


B\B, - B/B* 


B\B 9 - B/B* 


= I. 


dA 2 {t) 1 r dATit) 

[ dAf(t) \ " [ dAT*{t) J 


(99) 


According to 


we have 


—b' 2 b 1 + b; b .f -b' 2 b 2 + B{ by 

In our case that B is a 2 x 2 matrix, we have 

B\B 2 - B%B* = 0--BlB 1 +B\B* = 0. 

Moreover, we need to have the following relation: 

B\B x -B t 2 B* = -B\B 2 +bJb* = (Bl+B 2 ly )-(B 2 2x +B 2 2y ) = I. 

where B\ = B\ x + iBi y , B 2 = B 2x + iB- 2y . Thus, we can 
assume that 


■ dAnt) ■ 


y/k 

0 

a(t) 

dA°/'*{t) 


0 

y/ k 

a*(t) _ 


dA 2 (t) 

dA*{t) 


( 100 ) 


B lx + B i y = cosh(r) 2 ; B 2 X + B° 2y = sinh(r) 2 . 

Hence, we may write the matrix B in the following form: 

g _ cosh (r) cos (a) + zcosh(r)sin(o!) sinh(r)cos(/3) + zsinh(r)sin(/3) 

sinh(r)cos(/3) — zsinh(r)sin(/3) cosh(r)cos(Q;) — zcosh(r)sin(a) 












































































































Since £?,/ are 2x2 matrices, we have 


( rir 1 -/)- 1 


l 

2 — 2B\ X 


B\ - 1 

B* B* - 1 


Therefore, the last term on the right side of equation (103 1 
can be expressed as: 


S(B _1 - J) -1 = 


2 — 2cosh(r)cos(o:) 

, [" cosh (r) cos (ck) — 1 0 

M 0 cosh(r)cos(a) 


- 1 


2 cosh(r)sin(o:) sinh(r)cos(/3) + isinh(r)sin(/3) 

sinh(r)cos(/3) — isinh(r)sin(/3) —icosh(r)sin(a:) 

-f 0 
0 


2 — 2cosh(r)cos(o:) 
icosh(r)sin(o!) sinh(r)cos(/3) + isinh(r)sin(/3) 

sinh(r)cos(/3) — isinh(r)sin(/3) — 2 cosh(r)sin(o:) 


( 104 ) 


Substituting ( |104| > into ( |103[ >, we have 


daft) 


1 1 

2 

la#(t) 


1 —- 
2 J 


0 

dAi (t) 

0 

sfn 

dA*(t ) 


a(t) 

a#(t) 


dt 


2 — 2cosh(r)cos(a) 

zcosh(r)sin(a) 

sinh(r)cos(/3) — zsinh(r)sin(/J) 

a{t) 


a*{t) 


dt. 


sinh(r)cos(/3) + isinh(r)sin(/3) 
—(cosh(r)sin(a) 


(105) 


Therefore, the closed loop system wi th D PA and a static 
squeezer has the dynamical equatio n (|105| l. The difference 
between this closed loon system ( | 1 fe [i and the original 
uncertain quantu m syst em ( |98| ) is the addition of the last term 
on right side of (|105|i whicn corresponds to the controller 


Hamiltonian 11>. To realize our desired controller Hamilto¬ 
nian // ( as in (97 1 , the following relation is required: 


2 — 2cosh(r)cos(a) 

icosh(r)sin(a) 

sinh(r)cos(/3) — isinh(r)sin(/3) 


sinh(r)cos(/3) + isinh(r)sin(/3) 
—(cosh(r)sin(a) 


= —iJK = 


0 

-0.5 


-0.5 

0 


(106) 


When a = 0,/3 = 0,sinh(r) = — |,cosh(r) = |,k = |, 

the relationship ( |106[ ) holds. In this case, we may have B 

which satisfies the Bogoliubov 


matrix as B = 


condition B^ JB = J. 
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